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, Abstract. The two-neutron unbound ground state resonances of 26 O and 16 Be were populated 

^ ■ using one-proton knockout reactions from 27 F and 17 B beams. A coincidence measurement of 

3-body system (fragment + n + n) allowed for the decay energy of the unbound nuclei to 
be reconstructed. A low energy resonance, < 200 keV, was observed for the first time in the 
24 + n + n system and assigned to the ground state of 26 0. The 16 Be ground state resonance 
was observed at 1.35 MeV. The 3-body correlations of the 14 Be + n + n system were compared 
to simulations of a phase-space, sequential, and dineutron decay. The strong correlations in 
the n-n system from the experimental data could only be reproduced by the dineutron decay 
simulation providing the first evidence for a dineutron-like decay. 



1. Introduction 

The addition or removal of neutrons from stable isotopes has been shown in many cases to 
drastically alter the structure of the nucleus [1-3]. In moving towards the neutron dripline 
the typical shell gaps observed in stable nuclei have been shown to disappear and new magic 
numbers are observed. For example, in the dripline nucleus 24 a large energy difference was 
measured between the vlsin and vOd^^ orbitals indicating a new magic number at N = 16 [4, 5]. 
Determining how the shell structure changes for nuclei with extreme neutron-to-proton ratios 
is critical for providing stringent constraints on theoretical calculations. Such measurements 



have already been used to demonstrate the importance of three-nucleon forces in describing the 
properties of very neutron-rich nuclei [6]. 

Radioactive-ion beam facilities have made it possible to produce and study nuclei at and 
beyond the neutron dripline. The neutron dripline is defined experimentally only up to the 
oxygen isotopes, where 24 is known to be the last bound isotope [7-9]. The fluorine dripline is 
known to extend at least to 31 F but heavier isotopes have not been confirmed bound or unbound. 
A striking difference is observed in comparing the location of the neutron dripline for the oxygen 
and fluorine isotopes. While the 8 protons in oxygen are able to bind up to 16 neutrons, the 
addition of a single proton enables at least 22 neutrons to be bound. No published theoretical 
calculation has been able to simultaneously reproduce both the oxygen and fluorine driplines. 
Thus, measuring the properties of the unbound oxygen isotopes can constrain the theoretical 
calculations and help offer insight into the structure of the neutron dripline. Predictions of the 
two-neutron separation energy for 26 have shown large variations [6, 10-13]. An experimental 
measurement of the two-neutron separation energy would constrain the theoretical calculations 
and, thus, provide new insights into the changing structure of the neutron dripline. 

These exotic dripline nuclei can also exhibit unique types of decay which can provide 
additional insight into the structure of the nucleus. Goldanksy was the first to propose that 
the simultaneous emission of two protons could be observed given a scenario in which the 
intermediate state was positioned above the initial state [14]. Similarly, if the same conditions 
are found at the neutron dripline the direct emission of two neutrons would be expected. 
This scenario exists for the two-neutron unbound 16 Be, which should decay through the direct 
emission of two neutrons to 14 Be. The limiting cases for the simultaneous two-neutron emission 
would be a phase-space decay where the neutrons exhibit no correlation or the emission of 
a dineutron cluster. The more exotic dineutron decay is of great interest as multiple studies 
have shown indirect evidence for the presence of a dineutron in neutron-rich halo nuclei, such 
as 6 He and n Li [15-21]. However, a direct measurement of a dineutron-like emission has not 
yet been observed. Two-neutron unbound nuclei provide an excellent opportunity to observe a 
dineutron type decay since the emission occurs directly from the ground state in comparison to 
measurements of the two-neutron decay from bound nuclei where the system must be excited to 
overcome the two-neutron separation energy and, thus, the correlations may be altered relative 
to the ground state. 

In the following, we report on the first observations of the ground state resonances of 26 
and 16 Be. These results were first reported in Refs. [22] and [23] for 26 and 16 Be, respectively. 

2. Experiment 

The radioactive 27 F and 17 B beams were produced at the Coupled Cyclotron Facility at the 
National Superconducting Cyclotron Laboratory at Michigan State University. Two separate 
experiments were completed to measure the ground state resonances. In each case, a stable 
primary beam was fragmented on a thick Be production target producing a wide variety of 
isotopes. The A1900 fragment separator was used to select the desired secondary beam through 
the removal of the other reaction products and primary beam. The relevant information about 
each experiment is presented in Table 1. 

Measurement of the three-body decay (frag + n + n) required coincident detection of both 
neutrons and the residual fragment, which was accomplished using the modular neutron array 
(MoNA) [24, 25] and the large gap 4 Tm sweeper magnet [26]. The large area neutron detector, 
MoNA, is composed of 144 plastic scintillator bars. Light guides and photomultiplier tubes are 
attached to the ends of each bar for detection of the light produced from the interaction of the 
neutron(s). Each bar is 200 cm x 10 cm x 10 cm and are typically configured into 9 walls, each 
16 bars high. The time-of-flight (ToF) of the neutrons to MoNA is measured with respect to a 
scintillator placed in front of the target. The angle and energy of a neutron can be determined 



Table 1. Details of the experiments used to produce the unbound O and "Be nuclei. 



Unbound 
nucleus 


Primary 
beam 


Secondary 
beam 


Reaction 
target 


Secondary 
beam rate 


Ref. 


26Q 


140 MeV/u 
48 Ca 


82 MeV/u 
27 F 


705 mg/cm 2 
9 Be 


14 pps 


[22] 


16 Be 


120 MeV/u 
22 Ne 


53 MeV/u 
17 B 


470 mg/cm 2 
9 Be 


250 pps 


[23] 



from the interaction point in the scintillator bar and the ToF, respectively. The sweeper magnet 
bends the charged particles ~43° into a suite of position sensitive charged-particle detectors. 
Both the particle identification and kinematical properties of the charged fragments can be 
determined from the sweeper magnet detectors (see Ref. [27] for more details). 

From the measured energy and angle of the neutrons and charged particle, the invariant mass 
of 26 (M26 ) and 16 Be (Mi6 Be ) could be calculated. The 3-body decay energy could then be 
calculated as -Edecay = M26Q — M24 Q — 2M n in the 26 case, where M24 G is the rest mass of 24 
and M n is the neutron mass. In the two-neutron decays, the 2-body decay energy can also be 
examined and is calculated, for the 26 example, as -Edecay = M250 — Mi^q — M n , where M250 is 
the invariant mass calculated from the charged particle and the first detected neutron in MoNA. 

3. Monte Carlo Simulation 

A detailed Monte Carlo simulation has been developed to simulate the production and decay 
of the unbound nuclei. While a general discussion of the Monte Carlo simulation is provided in 
Ref. [28], the simulation of the neutron interactions in MoNA is of particular importance. In 
measuring the two-neutron decays, it is imperative that the simulation be able to accurately 
reproduce the multiple scattering of neutrons within the array. A single neutron can scatter 
twice producing a false 2n signal in comparison to two neutrons each interacting within the 
array providing a true 2n signal. The simulation must be able to correctly reproduce the ratio 
of false to true 2n events in order to allow for the resonance parameters of the experimental data 
to be extracted. 

The MoNA simulation is built upon the Geant4 framework [29,30], which allows for the 
tracking of each neutron interaction throughout the array. In order to validate the accuracy of 
the Monte Carlo simulation, a "clean" experimental data set is required for comparison. Ideally 
this would be a case where only a single neutron is impinging on MoNA. The measurement of the 
ground state of 16 B by Spyrou et al. [31] was chosen for this reason. A one-proton knockout from 
C produced the unbound 16 B, which immediately decayed into 15 B + n. For each detected 
15 B fragment there should only be a single neutron in coincidence. Therefore, any multiple 
interaction events observed in MoNA must correspond to the multiple scattering of the single 
neutron and should allow for a detailed comparison with the simulation. 

In Fig. 1 the multiplicity and energy deposited distributions from hits in coincidence with 15 B 
are shown from the experiment in comparison to the Geant4 simulation. The standard physics 
classes, referred to as G4-Physics, included in the Geant4 package are unable to reproduce 
either of the distributions. This was found to be related to the treatment of the inelastic neutron- 
carbon reaction channels [28]. In the G4-Physics simulation, the Geant4 cascade model is used 
to determine how the different possible inelastic channels are partitioned [32]. In comparison 
to the G4-Physics, a custom neutron interaction model, referred to as menate_r [33,34], was 
incorporated into the Geant4 framework. While the neutron(s) was still transported using the 
typical Geant4 methods, the menate_r class determined how each interaction was treated. 
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Figure 1. (Color online) The experimental (a) multiplicity and (b) deposited energy 
distributions are compared to the Monte Carlo simulation using the menate_r and G4- 
Physics classes. The simulated multiplicity distributions were scaled to match the number 
of multiplicity = 1 experimental events. The simulated deposited energy distributions were 
normalized to the total area of the experimental distribution. 



The menate_r code determines the partitioning of each of the different inelastic neutron-carbon 
reaction channels based on a compilation of experimental cross sections. Thus, the relative 
rates for each reaction channel should be accurately reproduced. A drastic improvement in the 
accuracy of the simulation is observed with the inclusion of men ate_r within Geant4 [28] , as 
shown in Fig. 1. 

As mentioned above, the two-neutron measurements are dependent on the ability to separate 
true and false 2n events. In both the 26 and 16 Be studies, causality cuts placed on the relative 
distance and velocity between the first and second hits in MoNA were used to remove false 2n 
events. In Fig. 2 the relative distance (D12), angle (#12), and velocity (V12) are shown from the 
experimental 16 B decay data [31] and the simulation using menate_r. Since each event in the 
experiment should only have one neutron in coincidence, the distributions, which require two 
hits in MoNA, are representative of multiple scattering or false 2n events. It is clear that false 
2n signals are largely correlated with a small -D12, roughly less than 30 cm, and a V12 less than 
about 10 cm/ns (which is equivalent to the beam velocity). This demonstrates the basis for the 
applied 2n casuality cuts used in the analysis of the 26 and 16 Be. The most important aspect is 
the excellent agreement shown between the experimental data and the Monte Carlo simulation. 
This indicates that the simulation is able to accurately reproduce the multiple scatterings of the 
neutrons within MoNA and should, therefore, be able to reproduce the effects of applying the 
2n causality cuts. 

4. 2e O 

As mentioned above, the unbound ground state of 26 was searched for using a one-proton 
knockout from a F beam and measuring the triple coincidence of the charged particle and 
two neutrons. While theoretical predictions have varied significantly on the prediction of the 
ground state position, the continuum shell model (CSM) has shown excellent reproduction of 
the measured properties of the oxygen isotopes and has predicted the 26 ground state to be 
unbound by only 21 keV [13,35]. Thus, the presence of low energy neutrons would be an 
immediate indication for the 26 ground state resonance. 

In Fig. 3(a) the two-body, 24 + n, decay energy is presented. A previous measurement of 
the 25 ground state, populated using a 26 F beam, reported a resonance at 770±fg keV [4]. As 
shown, the experimental Ed eca y spectrum does show an increase around 700 keV, which is likely 
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Figure 2. (Color online) Relative (a) distance, (b) angle, and (c) velocity between the first and 
second interaction in MoNA from the experiment is compared with the Monte Carlo simulation 
with the menate_r neutron interaction model. The simulated distributions were normalized to 
the total area of the experimental distributions. 



associated with the population of the O ground state. However, more interesting is the strong 
increase in E^ ecay below 400 keV, which is the first indication that the 26 O resonance has been 
observed. Fig. 3(b) shows the 3-body decay energy, where the 26 is reconstructed from the 
charged particle and the first two hits in MoNA without causality cuts applied. Again, a clear 
increase in the decay energy is observed below 500 keV indicating the presence of a low-energy 
resonance in the 26 O system. 

The Monte Carlo simulation was used to simultaneously fit Figs. 3(a) and 3(b). Three decay 
channels were considered in the Monte Carlo simulation as shown in Fig. 4. The one-proton 
knockout populating the ground state resonance of 26 (long-dashed red line) was modeled using 
a Breit-Wigner lineshape and assuming the two neutrons were emitted simultaneously with a 
phase space distribution. The population of the first excited state (short-dashed green line) was 
also considered and simulated as a sequential decay through the ground state resonance of 25 0. 
Lastly, the 25 ground state (dot-dashed blue line) could also be directly populated through a 
pn-removal from the 27 F beam. 

In the fit of the experimental data, the amplitudes of each decay channel were allowed to vary 
freely. However, since the data are not sensitive to the parameters of the possible contributions 
to the high-energy continuum the resonance parameters of the sequential decay from the first 
excited state and the direct population of the 25 O ground state were kept constant. The 25 O 
ground state resonance was simulated with an energy-dependent Breit Wigner with E = 770 keV, 
r = 172 keV, and L = 2 following the results of Ref. [4]. The first excited state of 26 was also 
represented by a Breit-Wigner with E = 2.0 MeV, taken from the CSM calculations [13,35], 
and r = 200 keV. The 26 ground state resonance was modeled as a L = 2 energy-dependent 
Breit-Wigner and both the resonance energy and width were free parameters in the fit. 

As shown in the top panels of Fig. 3, the sum of the three components in the Monte Carlo 
simulation (solid black line) were able to reproduce the experimental decay energy spectra. 
Through minimizing the x 2 °f the fit the 26 O ground state resonance energy was determined to 
be at E = 150±fg keV and was insensitive to the width of the resonance. The 26 O ground state 
resonance (red long-dashed line) is shown in Fig. 3 with E = 150 keV and T = 5 keV. 

The grey shaded region in the 3-body decay energy (Fig. 3(b)) represents the false 2n 
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Figure 3. (a) Decay energy spectrum of the two-body system, 24 + n. (b) Three-body 
decay energy spectrum for the 26 system, (c) and (d) show the three-body decay energy 
spectrum after the subtraction of the simulated false 2n events and after applying the causality 
cuts, respectively. The experimental data are shown as solid black circles. The three decay 
channels used in the Monte Carlo simulation are described in the text and shown in Fig. 4. The 
grey-shaded region in panels (b) and (d) represents the simulated false 2n component in the 
spectra. 
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Figure 4. Energy and decay level diagram for 26 decaying to 24 0. The three decay channels 
used in the Monte Carlo simulation are shown. The position of the 25 O ground state was taken 
from the experimental measurement of Hoffman et al. [4] and the 26 O energy levels are from the 
CSM calculations [13,35]. 
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Figure 5. Energy and decay level diagram for 16 Be. The solid blue lines represent previous 
experimental observations of the first excited state in 14 Be [37] and a lower limit on the position 
of the 15 Be ground state [38]. The 16 Be ground state measurement from the current work is 
shown as the solid red line. The dashed black lines represent shell model calculations using the 
WBP interaction [36]. 



component estimated by the Monte Carlo simulation. The false 2n events consist of a significant 
portion of the spectrum, yet the 26 ground state resonance appears to consist of mostly true 
2n events. Removal or reduction of the false 2n events can provide additional evidence for the 
presence of the ground state resonance. Two methods were used for this: (1) the simulated 
false 2n events were subtracted from the experimental 3-body spectrum as shown in Fig. 3(c) 
and (2) the causality cuts were applied to the data as shown in Fig. 3(d). After subtracting 
the simulated false 2n events, the spectrum contains only the ground state resonance and a 
background level of counts at higher decay energies. Applying causality cuts of Dyi > 25 cm 
and V12 > 7 cm/ns to both the simulated and experimental data shows a reduction in the false 
2n events by a factor of 3 with the ground state resonance observed well above the false 2n 
component. The persistence of the low decay energy peak after the removal and reduction of 
the false 2n events demonstrates strong evidence for the observation of the 26 O ground state 
resonance. 

5. 16 Be 

Shell model calculations, as shown in Fig. 5, have predicted that 16 Be is bound with respect to 
one- neutron decay and unbound with respect to two-neutron decay [36] . This provides a unique 
scenario in which the 16 Be should decay through the simultaneous emission of two neutrons. 
Experimental measurements (blue solid line) of the first excited state in 14 Be [37] and a lower- 
limit on the 15 Be ground state [38] both indicate that the shell model calculations appear to be 
reasonable, which suggests that the 16 Be ground state energy should be around 1.0 MeV. 

Based on the scenario presented in Fig. 5, three modes of decay for the 16 Be can be envisioned. 
The two-neutron emission could proceed via sequential decay through the tail of the broad l/2 + 
s-state around 3.5 MeV predicted by the shell model calculations. The two- neutrons could also 
be emitted simultaneously in an uncorrelated phase-space decay or as a correlated pair in a 
dineutron-like decay. Examination of the correlations in the 3-body decay can provide insight 
into the different decay modes. 

The 16 Be ground state was observed through the coincidence measurement of the 14 Be + n + n 
system. The 3-body decay energy spectrum is shown in Fig. 6(a) with causality cuts of 
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Figure 6. The three-body decay energy spectrum for 16 Be is presented in panel (a). Panels 
(b) and (c) show the two-body decay energy spectra in which the first and second time-ordered 
neutron was used, respectively. Panel (d) shows the decay energy of the n-n system, panel (e) 
shows the cosine of the opening angle between the two neutrons, and panel (f) shows the angle 
between the first neutron and the 14 Be fragment. The causality cuts discussed in the text have 
been applied to all the spectra, panels (a)-(f). 



D12 > 50 cm and V\2 > 11 cm/ns applied. A clear peak is observed in the 3-body decay 
energy spectrum around 1.5 MeV indicating the presence of the 16 Be ground state resonance. 
Additionally, the 2-body decay energy spectra are shown in panels (b) and (c) using the first 
(nl) and second (n2) time ordered neutrons measured in MoNA. 

The three different decay modes of 16 Be were simulated within the Monte Carlo simulation 
for comparison with the experimental data. In each case a Breit Wigner lineshape was used to 
represent the 16 Be ground state resonance. The correlations and partitioning of the energy 
between the emitted neutrons were then determined following the decay calculation. The 
sequential decay was calculated following the formulism presented in Ref. [13] for the emission 
of two / = neutrons with an intermediate s-state at about 3.5 MeV. The simultaneous phase- 
space emission was simulated using the calculation of Ref. [39] . In the dineutron simulation, 
the two neutrons are emitted together as a cluster and then immediately break-up. The break- 
up energy of the n-n system was calculated using a s-wave line shape [40] based on the free 
space n-n scattering length of -18.7 fm [41] and a pairing energy of 0.87 MeV [42]. The kinetic 
emission energy of the dineutron cluster is the difference in the total energy of the 3-body system 
determined from the ground state Breit Wigner and the break-up energy of the 2n cluster. 

In panels (a), (b), and (c) of Fig. 6 the best fit results of the Monte Carlo simulation for 
each decay mechanism are shown in comparison to the experimental decay energy spectra. As 



expected, the decay mechanism is insensitive to the shape of the decay energy spectra and the 
experimental data are fit well with each decay type. A x 2 minimization, varying E and V of the 
Breit Wigner, was completed for each decay mechanism. The resulting best fit corresponds to 
a "Be ground state resonance at E = 1.35(10) MeV with V = 0.8±g;^ MeV. As shown in Fig. 5 
(red line), the measured 16 Be ground state is in agreement with the shell model prediction and 
presents a scenario for the observation of genuine 3-body process. 

The correlations in the 14 Be + n + n system are shown in the bottom panels of Fig. 6. Panel 
(d) shows the decay energy of the n-n system. The opening angle between the two neutrons, 
in the frame of the 3-body system, is shown in panel (e), while the angle between the first 
neutron and the 14 Be is shown in panel (f). Unlike the E^ eca y spectra shown in panels (a)-(c), 
a clear sensitivity to the decay mechanism is observed in the correlations presented in panels 
(d)-(f). The sequential and phase-space decay mechanisms both show a broad E n+n distribution 
and a relatively flat cos(9 n - n ) distribution, which do not reproduce the observed correlations 
from the experimental data. In comparison, the dineutron decay simulation is able to reproduce 
the general features of the experimental data showing a peak in the E n+n distribution and 
a preference for small neutron-neutron opening angles. Overall, the agreement between the 
dineutron decay and experimental data implies that a strong dineutron-like correlation exists 
between the two neutrons emitted in the decay of 16 Be. While the decay is most likely not a 
true two-step process (emission of dineutron followed by the dineutron break-up), the presence 
of the n-n virtual state (a s = —18.7 fm) is clearly observed in this genuine 3-body decay process. 

6. Conclusions 

The ground state resonances of the two-neutron unbound 26 and 16 Be have been observed for 
the first time. The large area neutron detector MoNA and 4 Tm sweeper magnet were used to 
measure the charged fragment and two neutrons in coincidence allowing for the 3-body decay 
energy to be reconstructed. The 26 O ground state resonance was observed at E = 150±fg keV, 
in agreement with the sophisticated CSM predictions [13,35]. This extends the experimental 
measurements of the oxygen dripline and has already proven valuable in constraining theoretical 
calculations [43] , which have recently shown the need for three-nucleon forces to reproduce the 
measured 26 ground state. 

The 16 Be ground state resonance was measured at E = 1.35(10) MeV with a width of 
0.8^q' 2 MeV. This placed the 16 Be ground state below the experimental limit on the 15 Be 
ground state, which produces a scenario for a genuine three-body decay process. Three decay 
processes where simulated: (1) sequential decay through the tail of a broad s-state, (2) phase- 
space decay, and (3) emission of dineutron which proceeds to break-up based on the n-n free space 
scattering length. Only the dineutron decay simulation was able to reproduce the experimental 
3-body correlation observables suggesting that 16 Be decays through the emission of two strongly 
correlated neutrons. Since the simulated decay mechanisms represent extreme limits (dineutron 
versus phase-space), calculations using a more sophisticated 3-body model will be pursued in 
the future. Additionally, these results open the exciting prospect of finding other nuclei that 
exhibit similar signatures for a dineutron-like decay. 
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